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The axisymmetric collapse of a cylindrical air cavity in water follows a universal power law 
with logarithmic corrections. Nonetheless, it has been suggested that the introduction of 
a small azimuthal disturbance induces a long term memory effect, reflecting in oscillations 
which are no longer universal but remember the initial condition. In this work, we cre- 
ate non-axisymmetric air cavities by driving a metal disc through an initially-quiescent 
water surface and observe their subsequent gravity-induced collapse. The cavities are 
characterized by azimuthal harmonic disturbances with a single mode number m and 
amplitude a m . For small initial distortion amplitude (1 or 2 % of the mean disc radius), 
the cavity walls oscillate linearly during collapse, with nearly constant amplitude and 
increasing frequency. As the amplitude is increased, higher harmonics are triggered in 
the oscillations and we observe more complex pinch-off modes. For small amplitude dis- 
turbances we compare our experimental results with the model for the amplitude of the 



oscillations by Schmidt et al (2009) and the model for the collapse of an axisymmetric 



impact-created cavity previously proposed by Bergmann et al (20096). By combining 
these two models we can reconstruct the three-dimensional shape of the cavity at any 
time before pinch-off. 



1. Introduction 

The pinch-off of an axisymmetric air cavity in water is characterized by a finite-time 
singularity. The kinetic energy of the flow is focused into a vanishing small volume with 
a velocity whose magnitude diverges as the pinch-off moment is approached. Several 
experimental and theoretical scenarios have been recently considered in the study of this 



problem: a bubble rising from a capillary (Longuet-Higgins, Kerman fc Lunde|1991 



Oguz 



fc Prosperetti|[l993l |Burton, Waldrep fc Tab or ek| [20051 |Thoroddsen, Etoh fc Takehara 
2007 ), bubbl es in a co-flowing liquid QGordillo, Sevilla, R odriguez-Rodriguez & Martmez- 



Bazan 



2005 Bergmann, Andersen, van der Meer & Bohr 2009a), an initially necked 



bubble (Eggers, Fontelos, Leppinen & Snoeijer 2007), and cavities created through impact 
\ Bergmann, van der MeerT^ tijnman, Sandt ke, Prosperetti fc Lohse||2006| Gekle, van der 



Bos, Bergmann, van der Meer & Lohse 2008; Bergmann, van der Meer, Gekle, van der Bos 
fc Lohse|2009o| ). Depending on the case, the collapse might be initiated by surface tension, 
external flow, or hydrostatic pressure. However, irrespective of the cause, towards the end 
it is the inertia of the fluid that takes over in every case, and the collapse is accelerated 
as the radius of the cavity shrinks. 

The time it takes each of these systems to reach the inertial collapse regime varies 



by orders of magnitude (Gekle et al. 20096| ). Hence, it was not an easy task to deter- 
mine whether there was indeed a universal behaviour underlying this phenomenon. The 
first proposed model was a power law where the radius decreased proportionally to the 
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square root of the remaining time until collapse, r ( Longuet-Higgins et al. 1991 Oguz 



fc Prosperetti|[l993| ). Subsequent experimental and nu merical studies consistently found 



the behaviour deviated slightly from that \ power law ( Burton et a/.|2005 



Gordillo et al. 



20051 [Keim et a/.|[2QQ6l |Bergmann eTaI1[2006l (Thoroddsen et a/.||2007| |Bergmann eTal 



20096), generating doubts and starting a controversy about the universality of the phe- 



nomenon. Gordillo & Perez-Saborid (2006) and Eggers et al. (2007) theoretically showed 
how the power law varies weakly as a function of r due to a logarithmic correction. 

In conclusion, the axisymmetric problem converges to a universal self-similar solution. 
If axial symmetry is broken by a small azimuthal perturbation, a truly universal system 
would be expected to converge to the same solution . However, through experiments and 
simulations of an air bubble disconnecting from an underwate r nozzle, |Keim, M0ller~ 



Zhang fc Nagel|fl2006|); Schmidt, Keim, Zhang fc Nagel| p009| >; |Turitsyn, Lai fc Zhang 
(2009J); and Keim (2011) recently showed that a slight azimuthal asymmetry can trigger 
vibrations that persist in time. The fact that a small perturbation is not smoothed out 
indicates that the system possesses memory of its initial conditions. Here we conduct 
an experimental study of the evolution of azimuthal disturbances in the collapse of an 
otherwise axisymmetric cavity. 



Impact vs. detachment 

We study the cavity produced when a round disc with an azimuthal disturbance of 
its edge is driven downwards through the free surface of a water volume (see figures [I] 
and [2]). This is to be contrasted with the detachment of a bubble from a nozzle with a 
similar disturbance, which is initially determined by the competition of buoyancy forces 
with surface tension. The latter effectively smooths out large amplitude and high mode 
number (short-wavelength) perturbations, making nozzle experiments appropriate only 
for working with small amplitude, long wavelength disturbances. In addition, if the bubble 
is grown quasi-statically, viscosity can also play a role in smoothing perturbations due 
to the small Reynolds number of the water flow induced by the injected air. Therefore, 
this experiment does not allow much variation of the perturbations' mode number and 
amplitude. 

On the other hand, impact-created cavities are characterized by high Weber and 
Reynolds numbers from the beginning (provided that the collision speed is high enough) . 
Then, viscosity and surface tension play a marginal role in the formation of the cavity 
and are given no opportunity to erase features created by large-amplitude or high-mode 
perturbations. There is no initial surface-tension driven stage in the implosion; instead, 
the expansion of the cavity is opposed by hydrostatic pressure, which eventually starts 
the collapse, and is then quickly taken over by inertia. On top of this, since the cavities 
are created on a free surface, there is unobstructed optical access from the top, making it 
possible to track the shape of the horizontal section of the collapsing cavity. Disc impact 
experiments are consequently ideal for experimenting with the influence of geometry in 
cavity collapse. 

The effects of breaking the axial symmetry are clear (figure^). We present experimental 
results of cavities with disturbances of mode numbers 2 to 20 and amplitudes ranging 
from 1% to 25% of the mean disc radius. The experimental setup is described in ^2j 
We then explain the axisymmetric collapse model and the theory for the evolution of a 
perturbation in ^3] Our experimental observations are shown and discussed in Q The 
collapse of small-amplitude cavities viewed from the top and the side, and the comparison 
with theory are presented in §4.2| Finally, we discuss the collapse of high-amplitude 
shapes in §4.3| and draw general conclusions in ^5] 
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Figure 1: Collapse of a cavity created through impact of a round disc (a) and another 
using a disc with an azimuthal disturbance of mode number 20 (b) and amplitude of 
2% of the mean disc radius. In both cases the mean radius of the disc is 20 mm and the 
impact speed 1 m/s. The walls of the cavity in (b) acquire a structure that resembles 
the skin of a pineapple. This is explained by the oscillations triggered by the azimuthal 
disturbance. 



2. Experimental setup and procedure 

2.1. Setup 

The experimental setup consists of a linear motor that drives a vertical thin steel rod 
downwards through the bottom of a glass tank containing 50 x 50 x 50 cm 3 of water 
(figure [5| . A disc is attached horizontally to the top end of the rod with a mounting 
system that ensures that the disc is parallel to the undisturbed water surface and remains 
so throughout the impact. The motor is capable of a maximal 300 m/s 2 acceleration and 
its position can be controlled with a resolution of 5 |Ltm over its 1 m long track. Since the 
objects are not dropped into the water, velocity is a control parameter and not a response 
of the system. The setup allows for a precise control of the impact velocity on the range 



— 5 m/s. A more detailed description of the apparatus can be found in |Bergmann et al. 
( 20096] ). 



The shapes of the impacted discs are described by the functionr = S 'disci®) ? where 

Sdisc(0) = Rdisc + a m cos(ra<9), (2.1) 

Rdisc is the mean radius, a m the disturbance's amplitude and m its mode number (figure 
Discs are machined from a flat stainless-steel plate with a 2 mm thickness. The edges 
are sharpened to right angles with the intention of pinning the contact line to the lower 
edge and thus minimize the influence of wetting effects. Great care was taken to ensure 
that both the disc and the rod were dry before each run, as a single remaining drop on 
either is enough to noticeably alter the dynamics, especially towards the final instants 
before pinch-off. 

We recorded videos using a Photron SA1.1 high speed camera at frame rates from 
5,400 to 20,000 fps with resolutions ranging from 1024 x 1024 to 512 x 512 pixels. 
We image both top views of the collapsing cavity (focusing on the pinch-off plane) and 
side views to observe its structure. In order to avoid uneven optical reflections from the 
surface during top- view experiments, in some of the runs we diluted 1 g of milk powder 
per litre of water and shone light at the liquid, obtaining evenly scattered lighting. In 
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Figure 2: Experimental setup and examples of impacting discs. This camera- mirror con- 
figuration is used for top views; for side videos we simply move the camera down to the 
level where the light is. The shown discs have a disturbance amplitude of 10% of the 
mean radius (Rdisc)- 



this case the top surface of the disc was covered with black tape for improved contrast. 
For side views we used water without milk and diffuse illumination from the back. Top 
view videos were processed to extract the contour of the cavity in every frame and track 
the amplitude of disturbances as a function of time (j |4.2.1 ). Side views were directly 
compared to three-dimensional parametric plots of the modelling equations (j |4.2.3 ). 



2.2. Parameters 

The formation of a cavity after impact and its subsequent collapse are determined by 
the mean radius Rdisc of the disc, its shape, the impact velocity Vb, and fluid properties 
as kinematic viscocity z/, density p, and surface tension a. Hence, the dimensionless 
parameters of this experiment are the Reynolds number, Re = V^RdiscIv , the Froude 
number Fr = Vb 2 /(^disc#)? and the Weber number, We = pV^Rdiscj & - The surface 
tension (and thus We) is the only property that varies significantly following the addition 
of milk for some of the experiments. Since we focus on analysing the influence of the 
impactor's geometry, the impact speed and the mean disc radius were kept constant in 
all experiments reported here (Vb = 1 m/s and Rdisc — 0.02 m). In this way, the two 
variable control parameters -the mode number m and the amplitude OjfYi are related 
exclusively to the shape of the disc. In all of our experiments: Re ~ 2 • 10 4 , Fr ~ 5, 
and We ~ 300 — 400 (with measured surface tension values of [72 mN/m] for water and 
[47.1 mN/m] for the milk solution, respectively). These values indicate that the dynamics 
is dominated by inertia from the start of the experiment, and a scaling analysis reveals 
that this condition prevails throughout the experiment, making it unnecessary to consider 
their dynamic values. In the end, the only relevant control parameter during the evolution 
of the cavity is the Froude number since it is the parameter that determines when the 
cavity enters the inertial collapse regime ( [Gekle et a/]|20096[ ). 
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3. Models of cavity collapse 

3.1. Axisymmetric radial dynamics 
The model for an axisymmetric, impact-created cavity collapse by Bergmann et al. 



(20096) neglects vertical flow and its derivatives, assuming the collapse at each height z 
to be entirely decoupled from other heights, following previous works flLonguet -Higgins 
et a/.|[l99T Oguz & Prosperetti |l993). Since the flow is considered to be exclusively in 
the radial direction, the only relevant term from the continuity equation in cylindrical 
coordinates is: 



1 d 

r or 







(3.1) 



Integration of this equation, with the boundary condition that at the free surface the 
velocities of the water and the interface must be the same, i.e. u r (R) = R leads to: 



RR 

r 



which corresponds to a two-dimensional sink flow of strength Q(t) 
the radius of the cavity, and potential: 



(3.2) 

RR, where R(t) is 



$ = Q(t) ln(r) 

Conservation of momentum is expressed by Euler's equation of inviscid motion: 

1 dP du r du r 
p dr dt r dr 



(3.3) 



(3.4) 



which must be integrated from the cavity radius R to R^ . The upper integration limit is 
the length scale where the radial flow has decayed (R^ ^> R(t)); strictly, it should depend 
on the Froude number and time, but it is possible to determine an approximate constant 
average value from the experimental conditions and dimensions (Bergmann et a/.|2009fr| ). 
With AP being the (positive) difference between the hydrostatic pressure at R^ and the 
atmospheric pressure at the free surface of the cavity and accounting for the Laplace 



pressure jump (cr/R) across the interface at R(t), integration of equation (3.4) yields: 

"1 



AP 



o~ 
R 



-R z 



(3.5) 



where terms of 0(R/Roo) and smaller were neglected. 

In this way, the collapse at each height z is modeled like a two-dimensional Rayleigh- 
Plesset bubble collapse and the whole cavity is composed of a series of such collapses 



with different starting times as suggested by Lohse et al. (2004) in the context of the 
void collapse in quicksand. The original model does not include surface tension since an 
analysis of the dimensionless numbers from the problem reveals that surface tension never 
plays a major role in such a collapse. However, we have included it since disturbing the 
shape of the cavity creates regions of highly increased curvature where surface tension 
might play a role. 



During the inertial part of the collapse, the logarithmic term in (3.5 ) diverges as R goes 
to zero and thus the only way that equation can remain valid is by having the pre-factor 
of the logarithmic term go to zero. Integration from time t until the collapse time t co u 
yields a power law R(t) = a(t, 



coll 



(typical values found are 0.54 — 0.60) ( 


Burton et al. 2005 


Thoroddsen et a/.||2007 Keim 


et al. 


2006 


Gordillo et al. 


2005 Bergmann et al. 


2006 


20096) and theoretical studies 
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have shown that the exponent indeed has a weak dependence on the logarith m of the 
remaining collapse time, approximating to \ only asy mptotically at the end (Gordillo 
fc Perez-Saborid||2006| |Eggers et oT1|2007| [Gekle et a/.||2009frD. Nonetheless, the full 



theoretical result lies remarkably close to a power^law fit over many decades in time. 
Hence, we model the (dimensionless) mean radius R of our disturbed collapsing cavities 
as: 



R = a (t co u - t) 



(3.6) 



where the tilde indicates dimensionless quantities obtained by dividing length scales by 
R disc and time scales by Rdisc/V , i.e., R(t) = R(t)/R disc , t = tV /R disc and t co ii = 

tcollVo/Rdisc- 



3.2. Azimuthal disturbance 



Memory-encoding vibrations induced by a small geometric disturbance have been pre- 
dicted theoretically for any mode number m, and observed experimentally for cavities 
with m = 2 and m = 3 disturbances, namely bubbles released underwater from a slot- 



shaped nozzle, by Schmidt et al. (2009) and Keim (2011). The theoretical model was 



derived through a perturbation analysis of an azimuthal distortion to the geometry of a 



cavity with the behaviour described by equation 3.5 A brief explanation of the model 



follows for the sake of clarity. The complete derivation can be found in Schmidt (2008). 



Modelling the flow in an axisymmetric collapse as inviscid, irrotational, and incom- 
pressible (using Euler's equations) implies that there is no dissipation of energy. If so, 
the sum of kinetic and potential energies of the system will be conserved and can be 
expressed using the Hamiltonian: 



H(R,P R ) 



Pi 



2M(R) 



APirR 2 + a2irR. 



(3.7) 



The first term on the right is the total kinetic energy of the moving fluid, expressed in 
terms of the effective mass and its momentum, which are respectively: 



M(R) = 2 P ttR 2 ln(%) 
R 



P R = M(R) 



dR 
~dt' 



(3.8) 
(3.9) 



The second term is the potential energy due to the pressure difference between the 
fluid bulk and the cavity (at ambient pressure), and the third term is the energy cost 
of creating a free surface with the shape of the void. Applyi ng; H amilton's equations of 

We are thus faced with 



3.5 



motion R = OH/OPr and Pr = —dH/dR we recover equation « 
a dynamics with one degree of freedom, i?, and one constant of motion, namely the total 
energy. The implications of this are important: the problem is integrable, has a perfect 
memory, and according to the Kolmogorov-Arnold-Moser theorem if such a system is 
perturbed the new dynamics should closely follow that of the unperturbed situation; i.e., 
breaking the axial symmetry of the cavity by introducing a small disturbance of the shape 
should yield a collapse with the same leading order dynamics and new (approximately) 
conserved quantities, keeping it nearly integrable. 

Although in our experiments the disturbance of the shape of the cavity is characterized 
by a single mode m, the theoretical analysis considers a perturbation composed of a sum 
of Fourier modes cos(m#). For a conveniently chosen origin of the coordinate, the 
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perturbed shape of the void is 

S(0, t) = R(t) + a ™ CO cos(ra0), (3.10) 

m 

where R(t) is the mean cavity radius, which should follow the dynamics of the axisymmet- 
ric case, and a m (t) is the amplitude of each mode, which must be small when compared 
with the mean radius (a m (t)/R(t) <C 1) in order for the small perturbation theory to hold. 



Analysis of how the flow is modified by this shape disturbance ( |Schmidt]|2008| | Schmidt 
et g/.||2009 ) gives a linear second order ODE for the dependence of time evolution of the 



amplitudes a m on the mean radial dynamics: 

R , \ am(m 2 — 1) \ , 
i? (1 " m) + p& ) am = - ( 3 - n ) 

Equation (|3.11|) includes the additional influence of surface tension due to the surface 




disturbances on the right hand side of equation (3.10). We can find an approximate 
solution by substituting R(t) = a{t co u — t) 1 / 2 , neglecting surface tension, and solving the 
resulting Cauchy-Euler equation. In dimensionless form this gives 

a m (t) = a m (0) cos (\\/m - 1 In (t co u - t)) . (3.12) 

We see that the amplitude should oscillate with a constant magnitude and a frequency 
that diverges as t approaches t co u and that higher mode numbers will oscillate faster. 

Neglecting surface tension, a useful argument to understand the physics of the pre- 
dicted oscillations is the following: when the shape of the cavity is disturbed, its curvature 
is no longer uniform. Therefore, neither is the acceleration of the converging flow associ- 
ated to its collapse. As a result of continuity, convergence is stronger in regions with larger 
curvature, which consequently accelerate more and overtake the regions with smaller cur- 
vature (figures [3] and |4| , inverting the shape of the cavity. The higher the curvature, the 
quicker the overtaking becomes; thus with larger mode numbers more oscillation cycles 
are visible (figure [5| . 

Qualitatively, there is a connection to other instabilities that occur on accelerated fluid 
interfaces, like the Ray leigh- Taylor (RT) and Richtmeyer-Meshkov (RM) instabilities. 
Quantitatively this connection is less clear since, instead of being constant (RT) or shock- 
like (RM), in our case the acceleration is rapidly increasing in magnitude, and even 
diverges as r — » as a consequence of continuity. 

Equation ( |3.11 ), along with that for the axisymmetric radial dynamics (3.5) are the 



ingredients for the comparison of the observed cavity shapes obtained from experiments 



with theory (£4.2). 



4. Experimental observations 

4.1. Breaking the axial symmetry: general collapse mechanism 

Figure |3] shows the collapse of an elliptical cavity (which can be approximated as a m = 2 
disturbance to a circle) where the longer side of the cavity closes first. Initially, the shape 
of the cavity is the same as the impactor that created it (figure^). As it closes, the 
points that were originally farthest apart come towards each other at a higher speed than 
that of the end points of the minor axis (figure [3)3). Eventually it becomes clear that 
the shape of the cavity has inverted its phase with respect to the impactor (figure [3^) . 
This shift can be considered as an amplitude inversion of the original shape, described 
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Figure 3: (movie online) Top view of a collapsing cavity (m = 2, Raise — 20 mm a2 = 
0.2bRdisc) focused on the pinch-off plane. The cavity (neck) and the disc are initially in 
phase (a). Since the curvature is higher at the top and bottom of the neck, acceleration 
along the longer axis is larger. This changes the shape of the cavity as it collapses (b and 
c). Finally, the two points that were originally farthest apart come into contact first (d). 

as 5(0, t) = R(t) + a m cos(ra0) where R(t) is the mean radius and the amplitude, a, is 
originally positive and then becomes negative. 

For small amplitudes (1 or 2% of the mean radius), oscillations remain linear and with 
a nearly constant amplitude until very close to the pinch-off moment. However, the fact 
that we have a disturbance with constant amplitude in a shrinking geometry implies 
that the disturbance is actually growing with respect to the mean radius and hence it is 
bound to become of the same order of magnitude at some point. When this happens, the 
system starts developing higher harmonics and the linear oscillation model can no longer 
describe the events. At that moment, we say that the collapse evolves into a non-linear 
behaviour. This was observed for all mode numbers from 2 to 20. In §4.2| we will first 
discuss the linear regime, and afterwards turn to the non-linear effects in §4.3| 

4.2. Effect of small- amplitude disturbances 

4.2.1. Events on the pinch- off plane 

Small disturbances do not decay during the collapse of the cavities formed in our 
experiments. In figure [5] we show as an example eight snapshots of the evolution of an 
m = 16 cavity. Within the bounds set by the used frame-rate and resolution (20000 fps 
and a scale of approximately 86 /im/pixel), the cavities retain memory of the shape that 
created them throughout the whole collapse process. We tracked the edge of the cavity 
5(0, £) at every frame from the top- view videos, and found the mean radius R(t) and 
disturbance amplitude a m (t) at each time by fitting the curve: 

5(0, t) = R(t) + a m (t) cos(m0 + m (*)), (4.1) 
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(a) 

Figure 4: (movie online) Image of a m = 3, Rdisc = 20 mm, a = 0.1ir^ sc cavity, approx. 
10 ms after disc impact (a). Approximately 90 ms afterwards, it is clear that the shape 
has inverted with respect to the original (b). 



which is eq. (3.10) for a single mode perturbation. Since (f> m (i) was found to be nearly 
constant in time, we chose the reference angle = to be any line of symmetry in order 
to match with the model's description of the free surface. 

We determined the proportionality constant a and the power-law exponent f3 by fitting 



equation (3.6) to the experimental data of the mean radius. The values found for a were 

around 1.4 for all realizatio ns. F or a constant sink- flow of dimensionless strength Q = RR, 
with R given by equation (3.6) it follows that Q = /3a 2 r 2 ^ -1 , w here r = t co u — t. Using 
f3_ = \ (from the power-law model for the axisymmetric case 

2 



(1991]) and |Oguz fc Prosperetti| ( |1993[ )) we find that Q = ±a 2 . With a « 1.4 as found in 
our experiments Q exp « 1. The appropriate dimensional scaling of the flow in this work 
is Qtheo ^ RdiscVo = 0.02 (m 2 /s), which upon nondimensionalization with the scales 
Rdisc and Vb becomes Qtheo = 1, thereby confirming that the mean flow behaviour in 
our experiments is similar to the axisymmetric case. Values for f) were found to lie in the 
range 0.57 — 0.60, which is also consistent with previous works (see £3.1). 

Next we constructed the theoretical curves for the evolution of the amplitude a m by 
introducing the fitted power-law into equation (3.11). The initial condition for a m was 



determined from experimental observations and taken at a maximum of the curve; hence, 
the initial condition for the derivative is d m = 0. Figures [6] and [7] show the experimental 
results compared with the theory in semi-logarithmic plots of the mode amplitude versus 
the mean cavity radius, nondimensionalized using Rdisc ( m these plots time increases 
from right to left). The amplitude in the theoretical curves neither blows- up nor decays; 
it stays roughly constant in time. This is confirmed by the experimental data, at least 
during the first oscillations. The amplitude in experimental data drops at the end since 
we lose the capability to faithfully track the edge of the cavity towards the final collapse. 

The oscillation amplitude is roughly preserved, but as R(t) collapses, the relative dis- 
turbance a m (t)/R(t) grows (figure [8]). Since linear oscillations occur only for a(t) <C R(i), 
once this condition is not fulfilled non-linear effects overtake the dynamics, adding com- 
plexities to the shape of the cavity and increasing the difficulty of tracing its contour. 
The oscillation period looks constant, but as the horizontal axis is logarithmic, the fre- 



quency is actually increasing exponentially (chirping, cf. equation (3.12)). This is how 
the apparent contradiction between the universality of the axisymmetric system and the 
retention of the azimuthal disturbance at the same time, manifests itself in the dynam- 
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(a) r — 17 ms (b) r = 12 ms (c) r — 8.4 ms (d) r — 5.5 ms 




(e)r = 2.9 ms (f) r = 1.4 ms (g) r = 0.45 ms (h) r = 0.15 ms 



Figure 5: (movie online) Eight snapshots from a collapse with m = 16, Rdisc = 20 mm 
and a m = 0.02ii^ sc . The 16 peaks and valleys from the original shape can be clearly seen 
(a). The amplitude decreases on the way to inversion (b) and for an instant has a nearly 
round shape. The amplitude increases again (c) but it has now inverted with respect to 
(a). The process carries on (d,e) until we can still see a disturbance but cannot make out 
the details clearly at our experimental resolution (f,g). Finally, the void pinches off (h). 

ics. Part of the information (the mode amplitude) from the initial conditions is encoded 
and preserved -the cavity 'remembers' the shape that created it- but the chirping of the 
frequency makes it increasingly difficult to backtrack the evolution of the cavity as the 
collapse approaches, hence scrambling part of the information at the end ( |Schmidt et al. 
2009 ). 

4.2.2. The role of surface tension 

We shall now, as an intermezzo, comment on the role of surface tension. Accord- 
ing to the model, surface tension is expected to change the oscillation frequency and 
amplitude (figure [8]). Nonetheless, we see that experimental data acquired using water 
(°~w — 72 mN/m) and data with the milk solution (<r m = 47.1 mN/m) lay on top of each 
other (figures [6] and [7| . At first this might suggest that surface tension plays no role at 
all. If this were true, the theoretical curve without surface tension would fit these data; 
but it does not. Instead, the best fit is obtained with the surface tension of the milk 
solution. 

The milk powder dissolved in the experimental tank contains surfactant particles which 
adhere to the free and initially quiescent surface reducing surface tension. When the disc 
impacts and penetrates, a fresh free surface is rapidly created on the walls of the cav- 
ity, to which surfactants take time to adhere. Thus, for the duration of the experiment 
(~ 100 ms) surface tension must be effectively the same as water, or at least considerably 
higher than the value measured at the static surface. Nonetheless, from fit of the theoret- 
ical model we consistently find that the effective surface tension a e ff must be lower. We 
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Figure 6: Evolution of the dimensionless oscillation amplitude as the radius collapses 
(time goes from right to left). Results for several different mode numbers with an initial 
amplitude of 1% of the mean radius. Amplitude is expressed as a = a(t)/Rdi SC - Cir- 
cles: measurements with milk, crosses: plain water, solid line: theoretical prediction from 
equation |3.11| using the measured surface tension for the milk solution and for a given 
(dimensionless) mean radius R = a(t co u — i)P with a and f3 obtained from experiment. 



will now show that this can be interpreted as an effect of the axial curvature, neglected 
by the two-dimensional model. 

Considering the axial curvature as being related to the mean cavity radius R(t) (which 
is the curvature radius in the azimuthal direction) by a (nearly-constant) ratio 7 such 
that R ax (t) — — jR(t), we can write the Laplace pressure jump across the interface as 
o~ k = a (1/R(t) — l/jR(t)). The opposing signs in the curvatures are due to the hourglass 
shape of the cavity, which can be seen in figure [I] The right-hand side of the previous 
expression can be rewritten as a e ^/R(t) where a e ^ = a (1 — I/7), which for any 7 > 1 
gives a lower effective surface tension, thereby qualitatively explaining our experimental 
observations. Figure [l] also makes clear that R ax > R(t) and thus indeed 7 > 1. From 
side-view experimental images, we have observed that during the analysed part of the 
collapse, 7 varies between 2 and 4. The relation between the effective surface tension used 
for the theoretical fit (cr m ) and the value for water (a w ) is a m ~ |cr w , corresponding to 
7 ~ 3, which is in quantitative accordance with our estimation of the effect of the axial 
curvature. 
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Figure 7: Evolution of the dimensionless oscillation amplitude as the radius collapses. 
Results for several different mode numbers with a slightly higher initial amplitude of 
2% of the mean radius. Circles: measurements with milk, crosses: plain water, solid line: 
theoretical prediction. 

0.03^ ^ — — I ^ ^ I 




Figure 8: Theoretical behaviour of the amplitude, normalized by the disc radius, of a 
m = 20, <22o (0) = O.OlRdisc collapse with different surface tension values (left). — a = 0, 

Vm = 47.1 mN / m (milk), • • • a w = 72 mN / m (water). The same plot is shown on 

the right, but with the amplitude a2o(t) divided by R(t), to illustrate the growth of the 
relative disturbance. 
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4.2.3. The structure of the cavity 

We have so far tracked the instantaneous shape of the cavity on a horizontal plane at 



the pinch-off depth by looking from the top (see |Enriquez et al. 2010). Switching to a 
side-view allows us to see the complete evolution of the cavity at any given time before 
collapse in a single snapshot. Figure [9] shows four images of a cavity before it pinches 
off and one afterwards. We can see the walls of the cavity developing a structure that 
resembles the skin of a pineapple (see |Enriquez et al]|2011 ). Furthermore, the structure 
is not lost after collapse and is still clearly seen in the horizontal cavity ripples that form 



after pinch-off and were studied in detail by Grumstrup et al. (2007). 



We can reconstruct the shape by combining the model for axisymmetric collapse (3.5) 



with the equation for the perturbation's amplitude (3.11). First we solve the equation for 



the mean radial dynamics (3.5), and afterwards introduce the obtained R(t) and its time 



derivatives in (3.11). Since the models are two-dimensional and decoupled in the vertical 



direction, the three-dimensional shape is built by solving the equations simultaneously 



at several depths z as done in Lohse et al. ( 2004 ) for the cylindrical void collapse in dry 



quicksand. Figure 10 shows a parametric plot of the solutions just before pinch-off (10 i) 
and an experimental image at the same time (JToJd). We can improve on this result by 
using the same axisymmetric boundary integral code that was used in (Bergmann et al. 



2006 Gekle et al. 2008, 2009 a 



to obtain the undisturbed cavity profile R(z,t) which 
has been found to be in very good agreement with the experimental results (Bergmann 
et a/.||2009fr| ) and again use equation (13. 1 II) to superimpose the effect of the disturbance 
in exactly the same manner as described above. This procedure gives the shape in figure 
(flOfc) which is very similar to the experimental picture, capturing even small details. 



4.3. Effect of large- amplitude disturbances 

Increasing the amplitude of disturbances gives rise to more complex collapses. Such cases 
cannot be described with models derived from the analysis of a small perturbation, 
since there is a clear non-linear behaviour. Furthermore, the cavities' evolutions are 
not the same for all geometries. Phase-inversions are still observed and are explained 



by our continuity argument (£4.1) but it is very clear that the cavities no longer close 
approximately at a single point. Instead, a variety of closure types arise (cf. Turitsyn 
et al. 2009); for example: pointy and angular structures, finger-like forms, jets in the 
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Figure 10: (movie online) Side view images from experiments and model of m = 20, 
<22o = 2% collapses, (a) was obtained through the simultaneous solution of the model for 
two-dimensional axisymmetric collapse and the one for the amplitude of a disturbance; 
(b) is an experimental image; and (c) was obtained by applying the second model to a 
cavity profile that resulted from a boundary integral simulation of the axisymmetric case. 
This last method is more effective in reproducing the shape of the cavity, as it captures 
the axial curvature better. 



radial direction, and sub-cavities. Combinations of two or more of these events might 
take place. The symmetries of the cavities are always preserved, which means that the 
original single mode perturbation is dominant. However, other modes are excited and, 
since symmetry is conserved, the excited mode number m changes to an integer-multiple 
of its initial value. On top of this, there are situations in which one value of the angular 
position corresponds to multiple points of the cavity wall, thereby confirming that a 
simple mode description is no longer applicable. 

4.3.1. Some characteristic features 

In figure pT] we compare snapshots from m = 6 cavities with different initial disturbance 
amplitudes clq(0)/R(0) = 4%, 10%, and 25%. Each column is a time series of images 
belonging to one experimental realization. Corresponding snapshots in each series have 
been taken at equal remaining times to collapse (r). Whereas the 4% and 10% series 
are similar, with comparable features occurring at approximately equal times -albeit 
more pronounced for the larger initial disturbance- the 25% series differs considerably. 
A closer look at images aS and bS (r = 9 ms) in figure 11 reveals pointy, ridge-like jets 



being formed in singular cusps where the flow converges. For the largest disturbance, this 
happens only at r = 0.74 ms (c6). In time, such ridge-like formations can evolve into 
thicker, finger-like structures (b4)- 

We found these kind of shapes to be recurrent for most modes and amplitudes but the 
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time at which they happen varies in each case. When they occur too close to pinch-off 
there is no time for thickening and the collapse takes place with coalescence of the jets 
in the middle, which tends to happen when the initial amplitude is small (4% or less). 
Larger amplitudes seem to be more likely to form bigger sub-cavities. Still, the description 
cannot be generalized. For example, am = 6, a$ = 25% cavity will develop wide liquid 



structures that seem to crawl towards the center (figure 11 , cl) and eventually it assumes 
a shape with thin air arms (c4) through which water rushes in to invert the shape; in 
the end, the collapse (from the top) looks a lot like those with a small initial amplitude. 
In contrast, a m — 3, a 3 = 25% void will initially form similar thick structures, but they 



will come in contact with each other, forming large cylindrical sub-cavities (figure 12). 

The disturbance amplitude is roughly conserved during non-linear collapse. We have 
measured the longest and shortest diameters of cavities through several frames and from 
that estimated the mean radius and amplitude. The maximum amplitudes found have 



been of the same order as the impacting shape's. In figure 11 we can see that the os- 
cillation period changes noticeably when varying the initial amplitude. Look at the row 
where r = 9 ms and notice how the two cavities on the left (the smaller amplitude ones) 
have already started the shape inversion, as made clear by the jets that are being formed, 
while the rightmost void is still evolving towards the reversal. 

4.3.2. Sub-cavity formation 

By looking at the last images before the void closes it is easy to convince oneself 
that it is very unlikely that the cavities will pinch-off at a single point, thus sub-cavity 
formation must have happened in every case, although it might be extremely short-lived, 
not appreciable with our imaging capability, or both. The type of sub-cavity formation 
depends on whether the last inversion takes place early enough to allow for the thickening 
of the jet-like structures or not; and we have identified several different cases. For all 
mode numbers, m + 1 sub-cavities are formed -a central one surrounded by m satellites- 
except for m = 2, where there are just m of these (figure [3|. Sometimes the central sub- 
cavity is the smallest and collapses right away while the surrounding ones take longer to 
disappear. In other cases the exact opposite occurs. An example of the first case is shown 
in figure [12} where a m = 3, as = 25% void forms 4 sub-cavities. The center one is gone 
almost immediately but the other three live long enough to further partition and briefly 
become 6 little holes. When we can observe the sub-cavity evolution we always observe 
a phenomenon like in figure [l2b| where small jets are impinged from the contact points 
into all of the remaining sub-voids. 

4.3.3. Side view 

Another difference with small-amplitude collapses is revealed by looking at the pinch- 
off from the side. Figure [13] shows the pinch-off moment for a round disc and three 
m = 6 discs with amplitudes 4, 10, and 25%. The images in every case are just one frame 
before the cavity is definitely separated into top and bottom voids. Here the effects of the 
disconnection not taking place at a single point are manifested in drastic changes of the 
void's shape at the pinch-off point and accentuation of the top-down asymmetry, which 



affects the bottom part the most. Figure 13c shows three air columns with water around 



them; presumably there are three more columns hidden behind these (see figure [TT| b6), 
depicting a situation where seven sub-cavities were formed, the central one collapsed 



first, and the remaining six close afterwards simultaneously. In figure 13d we can see that 
the pinch-off took place in two stages; the "arms" of the cavity closed first, and smaller 
cavity is left in the center which collapses much like smaller amplitude voids (see third 



column in figure 11) 
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Figure 11: (movies online) Comparison of m = 6 collapses with different initial ampli- 
tudes. The oscillation periods of disturbances 4 and 10% of the mean radius seem a lot 
more similar to each other than to <2q = 25%, where it takes a considerably longer time for 
the shape to invert. Besides cusps (a6) and jets (b3,c6), we can see finger-like structures 
(b4, b5) that come in contact to form sub-cavities (b6). When the initial disturbance is 
very large, such features become very pronounced (series c). 
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(a) r w 12.8 ms (b) r w 5 ms (c) r « 3.5 ms (d) r w 0.5 ms 



Figure 12: (movie online) Example of sub-cavity formation for an m = 3, as = 25% 
collapse. Three large sub-cavities are formed around a smaller central one which quickly 
disappears. The inward pointing jets in (a) are due to the growth of the initial m = 3 
disturbance. The jets in (c) are due to the collision (shown in image b) and subsequent 
convergence of neighbouring regions of fluid which lead to six small cavities before closing 
(d). 




(c) m = 6, a 6 = 10% (d) m = 6, a 6 = 25% 



Figure 13: (movie online) Effects of axial symmetry breakup on pinch-off. Perturbed 
cavities no longer collapse at the same point and, as could be expected, the effects become 
more pronounced with larger amplitudes. Up-down asymmetry is noticeably accentuated, 
affecting more the void below the deep seal point. 
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5. Conclusions 

We have experimentally shown that breaking the axial symmetry of an impact-created 
cavity leads to oscillations of its walls as it collapses. A small, single-mode disturbance 
with an amplitude of 1 or 2% of the mean radius gives rise to linear oscillations of approx- 
imately constant amplitude and increasing frequency. The linear behaviour is maintained 
until the radius shrinks to a size comparable to the disturbance; afterwards, higher os- 
cillation modes evolve and non-linearity sets in. The mean radius evolves in the same 
way as the axisymmetric, universal case, making this system unique in the sense that it 
combines universal behaviour in the radial direction with memory of initial conditions 
in the azimuthal direction. Using two-dimensional models for the mean radius and the 
disturbance's amplitude we can reproduce the three-dimensional shape of the cavities 
and understand its structure. 

Increasing the amplitude of disturbances induces non-linear behaviour earlier in the 
collapse. The structures revealed in these cases attest to the beauty underlying collapse 
phenomena in fluid dynamics. We have observed a variety of pinch-off types arise and 
vary with mode numbers and amplitudes. Cavities preserve their symmetries in all cases, 
but the development of higher harmonics in this system is beyond our mathematical 
modelling so far. 
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